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Abstract 

^ ' The Lorentz-invariant gluon correlation functions {F^{x)F^{x')) and {FF{x)FF{x')) , cor- 

I responding to scalar and pseudo-scalar glueballs, are calculated for Kogan's and Kovner's 

( ■ variational ansatz for the pure SU (N) Yang-Mills wavefunctional. 

fT^ ' One expects that only one dynamical mass scale should be present in QCD; the ansatz 

I generates the expected scale for both glueballs, as well as an additional scale for the scalar 

■ glueball. The additional mass scale must therefore vanish, or be close to the expected one. 
This is shown to constrain the nature of the phase transition in the Kogan-Kovner ansatz. 

O ; 

PhI 1 Introduction 

Strongly interacting non-Abelian gauge theories require non-perturbative methods. One such, the 
, apphcation of the variational approximation to the Schrodinger formulation, goes back to the days 

^ ' of Feynman Though such methods are approximate, they incorporate many of the desirable 

. ^ , features of non-Abelian gauge theories, in particular the gluon condensate, asymptotic freedom, 

■ instantons, dynamical mass generation and confinement 1^-0. Variational methods have also met 
^ I success applied to Abelian theories |8|-[l0|]. 

■ - - ' In the Schrodinger formulation of pure 5[/(A^) Yang-Mills theory, the wavefunctional 5'[A^(a;)] 

is an eigenstate of the Hamiltonian 

H^l d'x^[Efix) + Bfix)], (1) 

where ^f(x) = tsAfir)^ ^ti^) = b^3k{d3Al{x)~dkA';{x)+gr^'^A)Al) and physical states satisfy 
Gauss' law as a constraint: 

[d^Etix) - gr^^A\EmA] - 0. (2) 
Kogan and Kovner proposed the gauge invariant Gaussian ansatz 

nAlix)] = J DU e^p[-^J d^xd'y AY^ix)S'^'5,,G-\x - y)A^'{y)] , (3) 

where under a local gauge transformation U{x), 

Alix) ^ A^'^ix) = S'^\x)A\{x) + Kix), (4) 



* b.gripaiosl@physics.ox.ac.uk 



1 



and DU is the SU{N) group invariant measure. Here S"'' ^ ^tT{T''UW'U), \f = Hrir^Wd^U) 

and ^ form an N x N Hermitian representation of SU (N): ^] = if"'^"^ with normahsation 
tr(T''r^) = 2(5°^. The A-field propagator has a mass parameter M, with respect to which the 
Hamiltonian is minimised. The authors found that M is of the correct size to generate the 
observed gluon condensate ]TT|] . 

In sections || and |^ of this paper, the Lorentz invariant gluon correlators (F^(a:)F^(x')) and 
{FF(x)FF{x')) are calculated using the same ansatz, and their scahng properties at large sep- 
arations are analysed. Naively, one expects that the correlators will decay as exp— M|x — x'\, 
with mass scale M. Firstly though, let us recall the details of the ansatz and calculation of the 
correlation functions therein. 



2 Correlator Rubric 

To calculate correlation functions, consider the generating functional 

Z = J DA ^[A]-9[A]; (5) 

the correlator of operator O is then 

iO) = ^J DA ^Ap^A]. (6) 

For gauge invariant O, one of the gauge integrations is redundant. Since the gauge transformation 
is linear on the gauge fields, the gauge field path integral is Gaussian. Using matrix notation for 
integration, viz. 

{ABmx, z) = Efc,, j d^y At^{x,y)B%{y, z), (7) 
the generating functional is proportional to 

J DADU exp[~^{A + afM{A + a)]exp[-^\'^{G + SGS'^)-^)X]. (8) 



Here, matrices are defined as 

St^ix,y) = S^\x)S{x-y)S,,, 
Gt!^ = G{x-y)5'^%„ 
M = S'^G-^S + G-^ and 
a^M-^S^G-^\. (9) 

Note that S is orthogonal whereas M. and G are symmetric. 

To calculate a correlator, substitute the operator in (|^) to obtain the correlator of some string 
of A fields. Then do the yl-ficld path integral; this is a free field theory. Thus 

{A'i{x)A]{y)) = a'i{x)a]{y) + M-'l^{x, y), (10) 

and so on via Wick's theorem. 

This done, one is left with an integral over gauge transformations U. The stratagem here is 
to factor the gauge transformation U into a product UlUh with parts dependent on high or low 
momentum modes of magnitude k: Uh with k > M and Ul with k < M respectively. Both the 
measure and action factorise similarly, since the group measure is invariant and the high and low 
modes are orthogonal. The A-field propagator takes the form G~^{k) = y/k^ + AP. For k > M, 
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take G ^(k) ~ k and treat the theory perturbatively. With U = exp[ig(j)°'T°'], to 0{g^) the high 
mode action takes the form 

Tff [0] d'xd'y d,r{x) G~\x ~ d.^^y). (11) 

One again obtains a free field, with propagator 

{^^\x)<^\y))^26'^\dtd\G'\x-y)\^^^^-^^1i'^''H-\x~y). (12) 

For k < M, take G^^{k) ~ M such that G^^{x — y)\k<M is local; the matrices Ul vary only 
on scales greater than 1/M and can be considered constant in the action]^ The resulting leading 
order low mode action is 

Tl[U] = ^tr j <fx d,Ul{x)d,UL{x). (13) 

The integration over high modes amounts to a renormalization group transformation. The effect 
on the low modes is to replace g in (^3|) by the coupling at the scale Af . This sigma model is 
calculated in the mean field approximation: the action is replaced by 

^-^tr J Sx [d.Ul{x)d,UL{x)+a^{UlUL - 1)], (14) 

subject to the constraints 

{uIUl) = 1, (15) 

{oUl) = 0. (16) 

So the Umn{x) act as N'^ free charged scalar fields with propagator 

{Ul^iix)UL,nniy)) = "^SkuSimi-dA + a'')-^5{x - y). (17) 
The momentum split allows the further simplifications 

S'^\x) = Sl\x), 

M''\x,y)^25^'^G-Hx-y), 
Xt{x)^Xt^ix)+XtHix), 

af{x) = ^X'i^{x) + lx'iHSi^{x). (18) 
Before tackling specific correlators, let us note also the following points: 

• The non-perturbative physics is captured in the low mode action; elsewhere we work to 0{g). 

• A factor 5*^ in a correlator contributes {U\Ul) in the Wick expansion, and therefore a factor 
of g^, from (|l^). The only leading order terms are those in which the orthogonality of the 
5'^, can be used, such that they cancel. 

• Orthogonality of the 5"^ can be used whenever they are separated by a scale less than 1/M. 
Thus, wherever, appears, it is sufficient to assume orthogonality, and integrate up to 
distance scale 1/M. Contributions on larger scales are higher order in g. 

• The action is quadratic in the decoupled Xl,Xh] correlators with an odd number of Al or 
Xh vanish. 



-"^In the action, (7^ are multiplied by G ^{x) ~ for \x\ > 1/M 
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• A correlator with an odd number of a and ergo A will therefore vanish to leading order. 

• Since Xf^ is just —dicj)^, f^ijkdj^kH ^ 

• \l is pure gauge and Abelian magnetic terms like {^iikdj^kh)^ vanish to 0{g). 

• Thus, any correlator with two or more eda terms must vanish to 0{g). 

• In seeking the scaling properties of correlators, constants and contact terms may be dis- 
carded. 

3 Scalar gluon correlator 

The Lorentz invariant gluon correlators are {F'^{x)F'^{x')) and {FF{x)FF{x')), where 

F^ = F^^F"'"' = 2{Bf ~ Ef), (19) 
FF = ^e^^.^F^^F^^ = iBtE^. (20) 

Discarding a contact term, the scalar correlator is thus 

(F^{x)F^{x')) = A(B^{x)B^{x') + E^{x)E^{x') - 2B^{x)E^{x')). (21) 

To 0{g), the chromo- magnetic term is 

Bfix)B-^\x) = e,,keun^e,,k^\ln^[^JAl^lAl,^]A^^^[A:}, 

+ gr^^d.AlJ^.A^J^A^d.A:^ + gr^^d',A^A^A^d,AldiAl,l (22) 

where a prime on a field indicates that its parameters carry primes. The terms of 0((/) have five 
A- fields and these vanish. Using (p^, the only terms with less than two e9as contain A^s only. 
To wit, 

{Bf{x)Bt\^))^^^keu.n^,,k^aJvd-,^^^ 

+d^d^-M-\i {x, x')dfd\-M~'Zln- 
+dp[^M-\i,{x,x')dfd'^-M-'^^l,{x,x')]. (23) 

Using eijkeum = SjiSkm - SjmSki and M'^^{x,y) = 25°-^5i-jG~^{x - y) yields (restoring primes and 
discarding a constant) 



+ i(A^^ - l)9jaj,'G(x - x')dpj:G{x - x') 

+ i(A^2 _ i)dpj'G(x - x')d}df,G{x - x'). (24) 



How do these terms scale? Recall 



G(x) = / ^ G{k)e^^-\ (25) 



where G{k) ^ = \/W + M^. Choose spherical polars with fca-axis along x. Then 

G(x) = —3- / dk , :sinfcx = —5— / dk e'''\ (26) 
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The contour can be deformed in C to surround the branch cut going from iM to oo (c/. pp. 
27-28]) and for large x approaches 

■..ji=:l,.jl«.,mx,hJL)"v-. (27, 



Similarly, 



where Ki^2{Mx) are modified Bessel functions ll3| . 

Carrying out the differentiation in (|2j), one finds that 

where the expected mass scaling is indeed present. 

Consider now the chromo-electric correlator. Here things become more complicated. Doing 
the functional differentiation in one obtains 

{Ef{x)E-^'{x)) = 6{N' - l)[G-\x - x')r 

-AG~\x~x') J d^yd'y' G'^y' ^ x')G-\y - x){Any)Any')) 

+ J (fySy'(fizd^z' G-^{y' - x')G-^{z' - x')G-^{y ~ x)G-^{z - x) 

{4{y')4{^')A1{y)A'i{z)). (30) 

plus constants. Using ( p^ and (|l^), and the orthogonality properties of the Sl, this becomes 

\[N^^l)[G-\x~x')f 
-l-G-\x ~ x') f dW G-\y' - x')G-\y - x)[X'i^{y')X'^^{y) + A^^(y')A^^ (jy)] 



+ ^ / d'^yd^y'd^zd^z' G-\y' - x')G-^{z' - x')G-\y ~ x)G-^{z ~ x) 



16 _ 

[\liy)KU^)KLiy)KL{^) + \H{y)\H{^)KH{y)KH{^) + 4A?//(y)^:L'^(y)A:2(^)A,'^(y)5^^(2/)A^^(z)]. 

(31) 

The inverse relation J d^y G'~^(x — y)G{y — z) — 5^{x — z) has been used. 

The averaging over the gauge transformations U remains. The orthogonality of the fc-modes 
ensures the following property: in any integral of the form J d^x f{x)g{x), where g contains only 
high (low) mode contributions, only high (low) mode parts of / contribute to the integral. Then, 
since G~^{x — y)\k<M — M5{x — y) is local, the integrations in ( ^ involving Al are trivial. 

The integrations involving \h are a little more tricky. Consider, for example, the simplest one: 
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Using (111), this is 



].G-\x~x') I d^ydy G-\y' -x')G'\y-x){\lMKH{y))- (32) 



-{N^-l)G-\x-x') J d^yd^y' G-\y' ~ x')G~\y - x)dfdf H-\y - y'). (33) 
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Fourier decomposition of the propagators yields the simple result 



-{N' - l)G~\x - x')G-\x - x')U>M. (34) 
The quartic term is similar. Since 

{^'^{y^''^iz)r{y)r{z)) - 4[(iV2 _ ifH-\y' - z')H-\y - z) 

+ (7V2 _ l)H-\y' - y)H-\z' - z) + {N^ - l)H-\y' - z)H-\z' - y)] , (35) 
the terms dependent on a; ~ x' become (Fourier decomposing again) 

These integrals can be evaluated, but are not very illuminating, amounting to Fourier transforms 
of various simple powers of k. If the theory is regularised with a cut-off, these are just spurious 
oscillatory functions of M\x ~ x'\. One obtains the same generic oscillatory result if one tries to 
calculate G~^{x — x')\k^nj explicitly; what is important for our purposes is that the only mass 
scale present is M. 

The remaining term in in ( |3l| ) is quadratic. The high mode propagator contributes a factor 
S''^ and the orthogonality property of the Sl can then be used. 

The low mode expectation values remain in (|3l|). The simplest terms are quadratic: 

(A^i(x')A^i(x)). (37) 

The Fierz identity for SU{N): t°jT^i = 2{SuSjk ~ ^S.jSki) implies 

KLW^Li^) - -^tTiuiix')d,ULix')Ulix)d.,ULix)). (38) 

The Ul contain modes with k < M; for x ^ x' , the unitarity of the Ul cannot be invoked. Doing 
the Wick expansion, one sees that these term are 0{g^). 
The quartic low mode term is (using ([l7|)) 

^{}:-^{x)KL{x)\'^^{x)\idx)) 

/Lf4 

=-—Air{d,U{x')d,U\x'))tr{d,U{x)d,U\x))) 



^i^i' ik.(x-x') 



This low mode term is very interesting. Consider the integral 

d^k e*'=-^ 1 /"^^ „ ke'^'^y 



(39) 



fc<M (27r)3 F + ^2 47T^y J_,j '^^ fc2 + ^2 ' (40) 

where {y = \x — x'\). In the limit M oo this is a simple contour integral. Closing the contour 
in the upper half plane, and enclosing the pole a± k = ia one obtains 

^e-y. (41) 
Airy 

Moreover, one can evaluate the integral for finite M using the exponential integral function (c/. p^ , 
p. 228]); one finds precisely this term, plus complicated terms involving M. But this term takes 
the form of a mass term; a appears as an additional mass scale in gluon correlation functions. 
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So now there are two mass scales, M and cr, in marked contrast to what one expects. How do 
we get ourselves out of this quandary? One way is to choose cr ~ M . Then the mass scale a is 
present and yet impossible to detect. The other way out is to force a to be zero. 

It seems then, that if one makes a somewhat arbitrary choice of cr = or M to fit the expected 
behaviour, then things work out satisfyingly. But these choices look distinctly less arbitrary in the 
context of Kogan's and Kovner's model. They consider the action ( p^ as a statistical mechanical 
system with SU{N) x SU{N) symmetry at temperature g'^{M). Since QCD is asymptotically 
free, large M corresponds to low temperatures, and conversely. One expects that the system 
will undergo a phase transition from an ordered state with spontaneous symmetry breaking to a 
disordered state with the full symmetry restored at some critical temperature. One is then left 
to speculate as to the critical temperature (coupling) and the nature of the phase transition. In 
the context of mean field theory, a acts as a mass gap; for a second order transition it vanishes 
continuously at the critical point, and is otherwise non-zero. 

Kogan and Kovner argue that the system's internal energy will be minimised at, or close to, 
the critical point, with a taking its critical value, corresponding to the mass gap. We have argued 
that cr can only take the values zero or M . Thus, the case cr = corresponds to no mass gap, 
and a second order phase transition, as predicted by the mean field theory. The case cr ~ M 
corresponds to a mass gap as big as the ultraviolet cut-off in the low mode model: clearly this 
corresponds to a strongly first order phase transition. Intermediate behaviour {e.g. a weakly first 
order transition) has observable consequences for the gluon correlators. 

Collecting up the terms generated by that is (|^, and (|39|), the chromo-electric 

correlator becomes, for cr = 

A{Ef{x)E[^\x)) = 6{N' - l)[G-\x - x')r 

- 4(7V2 _ i)c~\x - x')G-\x - x')\k>M 

(Pk Sk' '-'^2 



8A^^ / +{N^ ~ 1) 

k<M Jk>M. 



r-^(k\r.-'^(h'\ ^^'^'^ i{k+k').{x-x' 

(42) 



The final part of {F {x)F {x')) is (up to a contact term) 

-8{B'{x)EHx')) = -8(e.,.e.™ajA^arA"„z^^z^^), (43) 

where terms of 0{g^) or with an odd number of A-fields have been discarded. Calculating as 
before, we find that this is a contact term. 

After a good deal of effort, we have arrived at an expression for the correlator {F'^{x)F'^{x')) 
for cr = 0, namely the sum of ( p9| ) and (^). The explicit form of the correlator is not especially 
illuminating. What is significant is that unless cr = or M , there is an additional mass scale in 
the correlator. 



4 Pseudo-scalar gluon correlator 

Let us proceed to calculate the scaling properties of the pseudo-scalar glueball correlator. In 
calculating the scalar correlator, the a mass scale came from the quartic term in A^. Such a term 
is not present in the pseudo-scalar correlator, suggesting that only the expected mass scale M will 
be present. 

Discarding contact terms, the correlator can be written unambiguously as \Q{BB' EE') . Doing 
the functional differentiation in (^ and retaining terms even in A of 0{g) gives 

16(e.,fcaj^A^;.,9/4-[5-'5,,,G-i(x -x')- j Sy'Sz G-\y' - x')G~\z - x) X^\y) {z)]) . 

(44) 
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The first term contains two eds only, so its only non-zero part has (AA) A4 ^, i.e. 

16(iV2 _ i)G-\x - x')dJdfGix - x'), (45) 

with the expected mass scale M. The second term is more complicated. The possible non-zero 
parts in the four-point A correlator in (Q) are 

{d^Ald^^A,^A','^{y)Anz)) = aja/A^-i«^;(x,x')Al-^:';f 

+ d^M-'tt:{x,y')d-^M-";:,':ix',z) 

+ 2dJM-'l^:ix,y')d';a[^{x)at{z) 

+ dp'^^M-'r^,{x,x')a,%y)af{z). (46) 



The first term yields — 8(7V^ — l)G ^{x — x')djdj G{x — x') and the second term turns out to 
be a contact term. 



The third and fourth terms on the RHS of (46) have both high and low mode contributions 



For the low mode parts, the integrations in ( |44| ) are trivial, they result in a factor {Xl{x)Xl{x')) , 
which is 0{g^). 

The only remaining contribution comes from the high modes in the last term. This is 
16(iV2 _ ^^^QjQ. _ / ^G-i(fc)e^'=-("-" ) - dpj G{x - x')G-'{x - x')|fc>Af]. 

Jk>M (.^""j 

(47) 

Thus, the only mass scale in the pseudo-scalar glueball correlator is M. 

5 Discussion 

In the above, an extensive analysis of the gluon correlators in the Kogan-Kovner ansatz has been 
made. It has been shown that the gauge and Lorentz invariant gluon correlators contain the 
expected dynamically generated mass scale M. Moreover, it has been shown that an additional 
mass scale a is present in the scalar glueball, and should be observable, unless it vanishes, or lies 
close to M. This constrains the phase transition in the cr-model part of the Kogan-Kovner ansatz 
to be of second order, or strongly first order, respectively. If this is the case, then the scalar and 
pseudo-scalar glueball masses are degenerate. 
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